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Both wire- wound solenoids and cylindrical magnets can be approximately modeled as ideal, az- 
imuthally symmetric solenoids. We present here an exact solution for the magnetic field of an ideal 
solenoid in an especially easy to use form. The field is expressed in terms of a single function that 
can be rapidly computed by means of a compact, highly efficient algorithm, which can be coded as 
an add-in function to a spreadsheet, making field calculations accessible even to introductory stu- 
dents. In computational work these expressions are not only accurate but also just as fast as most 
approximate expressions. We demonstrate their utility by numerically simulating the experiment 
of dropping a cylindrical magnet through a nonmagnetic conducting tube and then comparing the 
calculation with data obtained from experiments suitable for an undergraduate laboratory. 



I. INTRODUCTION. 



The following article has been accepted by the American Journal of Physics. After it is published, it will be found 
at http://scitation.aip.org/ajp/ 

Solenoids and cylindrical magnets are staples of introductory physics laboratory experiments and demonstrations. 
When it comes time to put theory to the test, simple models for these objects are needed. An idealized solenoid - 
a solenoid with strictly azimuthal current in a thin sheet wrapped around a right circular cylinder - can serve as a 
reasonable model of an actual wire- wound solenoid, and can serve as an even better model of a permanent cylindrical 
magnet, provided that its magnetization is sufficiently uniform. 

At the introductory level, the magnetic field of an ideal solenoid of finite length can only be computed exactly along 
the symmetry axis, where the field can be expressed in terms of elementary functions. At off-axis points, geometry 
makes things more difficult and introductory-level students generally have no tools for obtaining even approximate 
values for the field except at very large distances where the field resembles that of a point dipole. This surely makes 
it harder for students to develop familiarity and confidence in dealing with magnetic phenomena. 

It is well known that the field due to a circular current loop can be written in terms of elliptic integrals, so by 
treating the ideal solenoid as a stack of loops, its magnetic field can be obtained by a straightforward integration. 
Alternatively, the field may be derived by solving a boundary value problem with cylindrical symmetry. In either case, 
exact expressions for the field have been known for over a century. They can be expressed in different forms using 
various special functions such as elliptic integrals, Heuman's lambda function, various Bessel functions, hypergeometric 
functions, etc. [H HI H 11 [S] 

It is true that some of these expressions can be rather cumbersome looking, and there are certainly many situations 
in which an approximate expression may provide better physical insight, be easier and/or less error prone to use, or 
serve some higher pedagogical purpose. However, when the expression for the magnetic field is merely being used as 
part of a computation of another quantity, it would seem that the only serious question about using an exact solution 
would be its computational speed. 

However, we easily found nearly a dozen recent papers[Hl EllinillllinilSlIIllISIISIIZlIIS] (see section |IV A] that 
included computations involving the field of a cylindrical magnet but did not make use of exact solutions in any form. 
The reasons for this choice probably varied from one paper to the next, but several papers explicitly stated that the 
theoretical treatment of the magnetic field due to a solenoid of finite length was too complicated for their purposes. 

There are apt to occur many occasions in which calculations involving solenoid fields arise. As a convenient tool for 
such situations, we present an exact solution in a form that (1) is algebraically less complicated, (2) does not require 
any previous knowledge of special functions, and (3) comes with a numerical algorithm that is simple and efficient. 
The field is expressed in terms of a single function, a generalized complete elliptic integral. This function is completely 
defined by an integral whose form occurs naturally in problems involving cylindrical symmetry. Numerical values can 
be computed by means of an algorithm that can be easily coded on even a programmable calculator or employed as a 
user-defined function or macro in a spreadsheet. In this form, the finite-length, ideal solenoid model is as simple and 
fast to use as the point dipole model in computations. 

This paper will present these exact expressions for the magnetic field of an ideal solenoid and exact expressions for 
its self-inductance, and will then provide a brief illustration of the effectiveness of these expressions by simulating the 
experiment of dropping a cylindrical magnet through a nonmagnetic conducting tube and comparing the calculations 
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with the results of some simple experiments. 

II. THE GENERALIZED COMPLETE ELLIPTIC INTEGRAL. 

In computing fields with cylindrical symmetry, certain integrals occur in a natural way. They are special cases of 
a function defined by a generalized complete elliptic integral (eel): 
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Appendix |A] describes code for an extremely efficient numerical algorithm for calculating values for C. The example 
code is presented in a dialect of BASIC that can be directly used as a user-defined function in spreadsheet program 
or can easily be translated to other programming languages. This appendix also contains further information about 
C, including the relationship of C to other forms of elhptic integrals and links to code. 

III. MAGNETIC FIELD EXPRESSIONS 

Consider a cylinder of length 2b and radius a wrapped by an azimuthal sheet of current Itotah equivalent to a 
tightly wound solenoid with a number of turns per unit length n carrying a current /, i.e. Itotai = 2bnl. The 
magnetic moment /i of the solenoid is /i = 2bnl ■ ira^. It is well known that along the symmetry axis of a such a 
solenoid the field, using cylindrical coordinates (p, ip, z) with the origin at the center of the solenoid, takes the form 
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which reduces to = i^onl for an infinite solenoid. For the general case, (see Appendix |b] for an outline of the 
derivation) the magnetic field components are: 
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FIG. 1: Field lines of an ideal solenoid with a length that is five times its diameter. 



These compact forms involve only a single function, C . They compute quickly and accurately both inside and outside 
the solenoid and are mathematically well-behaved except on the edge of the current sheet at p = a and z — ±6. 

These expressions reveal that, if distances are measured in units of a, then the magnetic field lines of an ideal 
solenoid depend only upon the ratio of length to diameter, b/a. Fig.[l]shows the structure of field lines for a solenoid 
with b = ba using the line integral convolution method employed by Sundquist[^ and Belcher. [7 Using the above 
expressions for the field, this image was produced by their JAVA program in a matter of seconds. The behavior of 
field lines within the solenoid contrasts sharply with that of external field lines and indicates why a single approximate 
formula in terms of elementary functions has difficulty representing the field at both near and far distances. 



Limits 



Along the axis of the solenoid (p = 0), k± — ^ — 1 and C(l, 1, 1, 1) — 7r/2, so Eq. Q reduces to Eq. Q. 
As 6 — * with 26n/ — Itotai remaining finite, a solenoid becomes a current loop and the field expressions above 
with < 6 << a do approximate those of a current loop. In the limit 6 = they take the form: 
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Finally, at large distances from the solenoid (r >> a, b), the field reduces to that of a point dipole: 
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B. Inductance 

It is worth noting that it is also possible to derive an exact expression for the mutual inductance of two coaxial 
ideal solenoids. |20 [ [2l ] [22] The self inductance L of a solenoid can then be obtained as a special case. In our notation 
the self inductance can be expressed very compactly: 

L= ^-^lo{naf [v/a2 + 62C(fco,l,l,2fc2) - a] , (16) 

where 

ko = , ^ (17) 

IV. FALLING MAGNETS 
A. Previous Work 

Faraday's law is often dramatically demonstrated by dropping a small, highly magnetized, cylindrical permanent 
magnet (radius a, length 26, mass m, magnetic moment [x) into a vertical, non-magnetic tube of conductivity cr, 
relative permeability 1, length L, inside radius r and the wall thickness be w << r) . [TTJ [TSl HSJ HSl [Ml Hi] The small 
masses and large magnetic moments of rare earth magnets give them a startling "hang" time. After the initial surprise 
subsides, students begin to ask questions. "How does the time of fall depend on the diameter or conductivity of the 
tube?" "How does it depend on the length of the magnet?" 

This experiment or similar ones have been analyzed in several papers. Pelesko et al. [8J and Roy et al. [S] use 
dimensional analysis to show that in a thin-walled tube the speed of the falling magnet should be proportional to 
{rngr^)/{a^^w), which demonstrates the power of dimensional analysis, though it is unable to provide information 
about how the speed might depend upon the geometry of the magnet. Other papers employ Faraday's law, but differ 



in the ways they model the cylindrical magnet. Hahn et al. [TT] derive Eq. (28 1 for the force on a magnet oscillating 
in a tube but then evaluate this by treating the magnet as a point dipole. Knyazev et al. [13j treat only the point 
dipole case, but expand the analysis to include high speeds not attainable in demonstrations. Levin et al. [TS] initially 
note that the ideal solenoid's external field is exactly equivalent to that of two uniformly (magnetically) charged discs 
at the top and bottom of the magnet, but then use point monopoles instead of discs in their calculations. Iniguez 
et al. [13] model the interaction of the magnet and the tube by means of an elaborate equivalent resistor network 
and provide a sample calculation in which the magnet is treated as a point dipole. Calculations based on the dipole 
approximation do not predict experimental results with much accuracy when the magnet fits closely within the tube, 
especially when the length of the magnet is greater than its diameter. 



After deducing Eq. (29 1 below, MacLatchy et al. [T7] model the cylindrical magnet as a stack of several polygonal 
loops and then compute its field from the Biot-Savart law. Although slightly cumbersome, this approach does offer 
adjustable accuracy. Partovi and Morris [18 offer a comprehensive treatment of a cylindrical magnet moving at an 
arbitrary, nonrelativistic velocity in an infinite tube of arbitrary thickness and permeability. This is a boundary- value 
problem with cylindrical symmetry that they solve, expressing the drag force on the magnet in terms of integrals 
involving Bessel functions with complex arguments. Their results are exact, though restricted to a steady state 
situation. The integral expressions in their paper are indeed daunting to look at, but the authors provide sample 
Mathematica code for computing them. 



The expressions for magnetic field presented in section III are extremely well suited to this problem 



B. Theory 

Previously cited references [TTJ [17] [18] discuss the theory in some detail, so we only provide a quick sketch here. 
Choose cylindrical coordinates with the z-axis vertical and the origin located (momentarily) at the center of the 
magnet. In an experimental situation, the conductivity of the tube can be determined by measuring its resistance per 
unit length R^: 

Rl = -7^, (18) 

where f = r + w/2. 
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FIG. 2: Geometry for a magnet magnet falling though a non-magnetic, conducting tube. 



As the magnet falls, the changing magnetic field within the tube walls is accompanied by an electric field that 
drives currents that cause ohmic heating. Although currents in the tube in turn induce currents within the permanent 
magnet, it is easy to show that under our experimental conditions the only significant energy losses are those within 
the tube walls. The speed of fall is so slow that air resistance is also quite negligible. In the following discussion, 
we assume that the magnet fits closely enough within the tube walls that its axis remains vertical and cylindrical 
symmetry is maintained during the fall. (When a small diameter magnet falls within a much larger diameter tube, 
one observes the axis of the magnet precessing about the vertical during the fall.) 

The electric field within the tube can be deduced by arguing that, in the reference frame of the falling magnet, 
there is only a magnetostatic field, so in the frame of reference of the tube, where the magnet has velocity v, there 
must be an electric field E = — v x B, or E^p = —VzBp. Alternatively, the field can be obtained from Faraday's law 
by considering a horizontal circular loop of radius f lying within the tube walls at height z' above the center of the 
magnet. Defining the upward magnetic flux though such a loop as 

f 

$(z')- J Bz{p,z')27rpdp, (19) 



the emf around the loop is 

emf - - '^^ ^^''^ - w (20) 

dt dz' dt dz 



since 



dz' 



If you then take two such loops separated by a vertical distance dz, you can visualize them as the edges of a small 
cylindrical Gaussian pill box. Since there are no magnetic monopoles, the total magnetic flux leaving the closed 
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surface of the box must be zero: 

$ (z + rfz) - $ (z) + Bp (f, z) 2TTfdz = = d<^ + Bp (f , z) 2'Kfdz. (22) 



Combining this with Eq. ( 20 1 , we have as before 

emf ^ ~{2TTfV:,)Bp^ E^27Tf. (23) 

The force acting on the faUing magnet can be deduced from energy considerations. The electric field within the 
tube drives currents that dissipate energy at a rate per unit volume of aE'^, so the total power lost P can be obtained 
by integrating this over the volume of the walls. If the walls are thin, the power loss when the magnet is at height z 
above the bottom of the tube is 

P = J (TEl2T:wfdz = 2TTWfavl j B^ (f, z') dz' = -V^Fdrag- (24) 

Alternatively, the Lorentz forces acting on the currents can be calculated directly. A horizontal slice of tubing of 
height dz located at height z above the center of the magnet is a circuit with electrical resistance 

1 27rr 

Re^ ^, (25) 

a w dz 



so the current due to the emf around the ring is 



GTnf 

dl = — — = —aVzW Bpdz. (26) 
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The vertical force exerted on this current by B is: 

dFz = -27rf dl Bp = a 2nfw v^BI dz, (27) 
and the force on the magnet then follows from Newton's third law. The total vertical force on the magnet is the sum: 



Fdrag^-^ J Bl{f,z')dz'. (28) 



in agreement with Eq. ( 24 1 . 

The magnetic field computes so quickly using Eq. ^ for Bp that this drag force can be used to numerically integrate 
the equation of motion for the falling magnet even though an integration must be performed at each time step. 

The powerful but light-weight magnets used in demonstrations reach a constant velocity within a fraction of a 
second. Since the magnitude of Bp declines extremely rapidly with z, this terminal velocity can be calculated by 
equating Fdrag to mg and, with little error, setting the limits of integration to ±(X): 

—mgRL 

^terminal — "Foe „n /_ r~7 ■ l^^j 

J-oo Bp ^ ) dz' 



C. Experiment 

We checked these predictions experimentally in the simplest possible way using equipment and procedures available 
in a typical undergraduate laboratory. We used copper plumbing tube (Table which is about 99% pure copper. 
Its magnetic permeability was not measured but taken to be the same as pure copper. The resistance per unit length 
was determined by running a current of several amperes through the tube while measuring a few millivolts potential 
difference across it. The conductivity of pure copper is about 110% of the standard lACS value, 5.8108 x 10^ S/m at 
20 °C, while the conductivity of typical copper tubing for plumbing is typically about 85% of the lACS value. 

We obtained six cylindrical magnets with various lengths but with the same half-inch (12.7 mm) diameter. We 
determined the magnetic moment fi for each one by using a small magnetometer to measure the field strength at 
several points along its axis and then adjusting the value of (nl) in Eq. ^ to give a best fit. 

Finally, holding a magnet vertically by its upper 4 millimeters, we inserted it into the top of the tube and released 
it while manually starting a timer that stopped when the bottom millimeter of magnet activated a photogate placed 
at the bottom of the tube. The results are summarized in Table [llj In addition to the measured average velocity of 
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Physical 


Electrical 


L = 1.478 m 


length between clips = 1.475 m 


r = 7.25 mm 


I = 4.95 A 


w = 0.7 mm 


AV = 4.49ml/ 


mass = 434 g 


a = 56.0 X 10^ S/m 




Rl = 5.37 X 10""*J7/m 



TABLE I: Physical and electrical properties of tube. Conductivity was determined by attaching a power supply to the ends of 
the tube with alligator clips, running a current I through the tube, and measuring the potential drop AV. 



No. 


b/a 


m(g) 


lj,{A ■ m?) 


"^average ij^l ^) 


ut-Partovi 


Wt-Eq. (291 


"^average 


1 


1.0 


12.1 


1.76 


0.0687 


0.0669 


0.0670 


0.0674 


2 


1.5 


17.9 


2.36 


0.1045 


0.1050 


0.1050 


0.1058 


3 


2.0 


23.8 


3.23 


0.1275 


0.1243 


0.1243 


0.1254 


4 


3.0 


36.4 


5.00 


0.1825 


0.1710 


0.1711 


0.1731 


5 


4.0 


48.2 


6.37 


0.2473 


0.2451 


0.2451 


0.2486 


6 


1.0 


12.9 


1.17 


0.1513 


0.1615 


0.1616 


0.1622 



TABLE IL Six cylindrical magnets with half-inch diameters and various ratios of length to diameter {b/a) dropped from rest 
through a length of copper tube. Average speed of fall is compared to predictions of terminal velocity using formula of Partovi 
and Morris or Eq. (29 1. Last column is the predicted average velocity obtained by numerically integrating the equations of 



motion of a magnet experiencing the drag force given by Eq. ( 28 1 



each falling magnet, the table shows pairs of terminal velocities numerically calculated both from Eq. (29 1 and from 
the lengthy exact expression for the terminal velocity in Partovi and Morris. (Both calculations were carried out in 
Python using the SciPy libraries.) The numerical agreement between these two different ways of calculating terminal 
velocity validates the approximations used in deriving Eq. ( 29 1 . 

In addition to calculating terminal velocities, we used Eq. (28) for the drag force on a falling magnet to numerically 
solve the equations of motion for z{t) and thus the total time of fall. Because of their differing lengths, the center 
of each magnet started out at a different height in the tube resulting in slight differences in their initial behaviors, 
which are plotted in Fig. [3] which shows that in all cases terminal velocity was achieved within a few tenths of a 
second resulting in terminal speeds that are only slightly smaller than average speeds of fall. Since the experimental 
procedures were not very sophisticated, the excellent agreement between theory and experiment was surprising. The 
experiment is clearly within the capabilities of introductory students. How much of the numerical simulation would 
be appropriate to ask them to do? That depends upon the level of the class and the amount of time allotted to 
this project. It is possible, though, to perform the computation of the terminal velocity at an introductory level if 
students are comfortable with spreadsheet programs or elementary programming concepts in, for example, VPython. 
Three things permit this. First, eel can be added to a spreadsheet as a user-defined function (see Appendix [A| , so 
Bp can be computed and graphed. Second, Bp is a symmetrical function of z, so the integral in Vterminai only needs 
to be evaluated from to oo. And third, at large enough values of z, Bp approaches the dipole form (Eq. (14)), 
so Bp ~ allowing the integral to be truncated at relatively small values of z. In fact, students can plot BJ^in 
a spreadsheet and see that it starts at at z = 0, reaches a sharp maximum near z = b, and decreases by several 
orders of magnitude by z = 36. As a result, the entire integral may be computed as a simple Riemann sum from 
z = to about z — 3b using roughly a step size dz ~ 0.015 - for which a spreadsheet is an ideal tool. Any concern 
about truncation errors could be addressed by calculating the remainder of the integral analytically using the dipole 
approximation for Bp. More advanced students, of course, might be acquainted with other programs or numerical 
methods for handling these tasks, but the elementary method described gives results within about a percent of the 
more sophisticated tools and can easily be improved by using more steps and/or smaller size steps. 



V. COMMENTS 



We have presented very compact expressions for the magnetic field of an ideal solenoid and have demonstrated that 
in numerical work they are easy to use, remarkably fast and can readily be incorporated into spreadsheet calculations. 
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FIG. 3: Speeds of magnets dropped through a copper tube. Each magnet was inserted into the tube and then released from 
rest. These are numerical solutions of the equation of motion including the drag force during the first fifth of a second after 
release. 



This makes it possible to easily simulate a variety of situations involving cylindrical magnets. In the case of the 
falling magnet demonstration, we have shown that a simple treatment of the problem provides results that agree with 
those from a more complicated analysis and are consistent with simple measurements. It should be clear that the 
methods used here can be applied to other cylindrically symmetric situations such as the electric fields of uniformly 
charged rings or cylindrical shells. Such fields can then be written in terms of eel conferring upon calculations all the 
advantages that have been described here. 

We hope that these expressions for the magnetic field will help dispel the notion that exact expressions for the 
field of an ideal solenoid are necessarily complicated or slow to work with and we hope they will encourage student 
investigations of magnetic phenomena. 



APPENDIX A: GENERALIZED COMPLETE ELLIPTIC INTEGRAL 

The generalized complete elliptic integral in Eq. ([T]) can be efficiently computed using an algorithm by Bulirsch |28j 
based on work of Bartky [29] who extended ideas of Landen and Gauss. This algorithm converges so quickly that, 
unless kc << 1, only 3 or 4 passes are necessary. The code shown in Fig. |4]is in a version of BASIC that can be loaded 
as a user-defined function into the Calc spreadsheet program (part of the free Open Office [35] software suite) where 
it can then be used as a normal spreadsheet function, making C accessible to non-programmers. Similarly, Microsoft 
Exeel allows user defined function coded in Visual BASIC. This code is simple enough to be treated as pseudocode 
that can be easily adapted to other languages. Alternatively, code in FORTRAN or C can be found on the internet [26] 
or in the first edition of Numerical i?ecipes. [3U] Code in FORTRAN, C, or Python is also available from the authors. 
Later editions of Numerical Recipes no longer mention eel and adopt instead Carlson's more general approach, which 
applies to incomplete as well as complete elliptic integrals. In terms of Carlson's functions, Rp and Rc, eel is: 

C{k,,p,c,s)=cRF{{),kll)+^-{s-pc) Rj{0,kll,p). (Al) 

Code for computing Carlson's functions may be found in the current editions of Numerical Recipes and elsewhere. 
However, for complete elliptic integrals, Bulirsch's algorithm is more compact and also has some other advantages. |27j 
This generalized complete elliptic integral includes all three standard Legendre forms as special cases: 

if(fc) = C (fee, 1,1,1), E{k) = C{k„l,l,kl) , n(n,fc) = C (fee, n+ 1,1,1), (A2) 

where 



fc= ^1 -fc2. 



(A3) 
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REM macro for Calc spreadsheet 
FUNCTION cel(kc, p, c, s) 
IF (kc = 0) THEN 

eel = "NaN" 

EXIT FUNCTION 

ENDIF 

errtol = .000001 

k = ABS (kc) 

pp = p 

cc = c 

ss = s 

em = 1 . 

IF (p > 0) THEN 
pp = SQR{p) 
ss = s/pp 

ELSE 

f = kc*kc 
q = 1. - f 

g = 1. - pp 

f = f - pp 

q = q* (ss - c*pp) 

pp = SQR( f/g ) 

cc = (c - ss) /g 

ss = - q/ (g*g*pp) + cc*pp 

ENDIF 

f = cc 

cc = cc + ss/pp 
g = k/pp 

ss = 2* (ss + f*g) 
pp = g + pp 

g = em 

em = k + em 

kk = k 

WHILE ( abs ( g - k ) > g*errtol ) 



k = 2*SQR(kk) 

kk = k*em 

f = cc 
cc = cc + ss/pp 
g = kk/pp 
ss = 2* (ss + f *g) 

PP = g + PP 

g = em 

em = k + em 



eel = (pi/2.)*(ss + cc*em) / ( em* (em + pp) ) 
END FUNCTION 



WEND 



FIG. 4: Algorithm for the generalized complete elliptic integral C{kc,p, c, s) coded in a version of BASIC for use as an add-in 
function in the Calc spreadsheet. 



In algebraic work, the following identity is useful: 




(A4) 
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FIG. 5: Geometry of an ideal solenoid showing notation used in applying the Biot-Savart law. 



APPENDIX B: DERIVATION OF SOLENOID FIELD 



The magnetic field of an ideal solenoid can be computed directly from the Biot-Savart law. The necessary algebra is 
only slightly more complicated than that which is needed to derive Eq. ([2| , which is commonly presented to students 
in introductory courses. The notation used for a±, k±, etc. is the same as that of section |T] 

First, the surface of the solenoid is divided into circular strips of width dz' as in Fig.jH The current in such a strip 
is nidz' . To calculate the field at a point x, we apply the Biot-Savart law to this circular loop and then add up the 
fields of the stack of strips that form the solenoid surface: 

b 

-6 

where 

R = X - x', (B2) 
and points along the strip are specified by the position vector 

x' = acos lys' i -I- asin j -I- z' k, (B3) 
and an infinitesimal step taken along the strip is 

dn' = (^—as'mip' i + acosip' dip' . (B4) 

Because of the cylindrical symmetry, we are free to choose coordinates in which x lies in the x-z plane, causing 
By{x.) to integrate to and allowing i?x(x) to be identified with _Bp(x), suggesting the notation 

x^pi + zk. (B5) 

Since 



R = (p — a cos if') i — a sin ip'^ + (z ~ z') k, 



(B6) 
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and 

dx' X R = a d(f' (z — z') cos Lp'i + (z — z') sin lyj'j + (a — pcos (f') k 
the field can be written as 



B(x) 



, f fionIa\ f I (z — z') cos (ys'i + (a — pcos 1^9') k 
dz ( — 1 . 



I — 2ap cos Lp' + + {z — z')^^ 



V2 



Integration over z' is elementary: 

-B„a 



dtp' cos 



1 



1 



z\ + + a? — 2ap cos ip' J z'^ + p"^ + a? — 2ap cos ip' 



and 



(x) 



Boa 



dip' 



(a — pcos ip') 



(p-^ + — 2a/9COS p') 



-^^i^+p^+a^— 2ap cos c/j' 



Q ^ — 2ap cos cp' 

To put these expressions into a form resembling eel, introduce a change in integration variable: 

2ip = n — ip' 

and, after using some trigonometric identities, observe that 

z^ + + — 2apcos (p' — z^ + (p + a)^ (cos'^ tp + k\ sin^ . 
The radial component of B then becomes 

^/2 

Bp (x) = Bo j d^ (cos^ V - sin^ tp) 



^ cos^ tp + k'^ sin^ t/i ^ cos^ -p + k"^ sin^ t/i 



(B7) 



(B8) 



(B9) 

(BIO) 

(Bll) 
(B12) 

(B13) 



Upon comparing each of the two terms in this integrand to the integrand in the definition of eel, we recognize that 
the radial field can be identified as Eq. ([3|. 

Similarly, the longitudinal component of the field becomes 



B, (x) 



Boa 



V2 



di) 



cos^ + 7 sin ip 



13- 



(P+'')J cos2 + 7' sin^ ^ I ^^^g2 ^ + fc^ sin2 ^ y^cos2 7/^ + fc2 sin^ ^ ^ 



(B14) 



which can be recognized as Eq. d4| 



In order to determine the magnetic field of a current loop (0 < 6 << a), it is simplest to return to Eq. (B9l 



and Eq. (BIO I, and, treating 6 as a small quantity, expand the integrands to first order in b and then repeat the 



transformations in Eq. (Bill and Eq. (B12) to put the integrals in forms resembling cel. The results are given in 
Eq. (ITTl) and Eq. (fT2l). 
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